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1.  Foreword 

There  were  two  major  results  in  this  project.  First,  we  elucidated  the  kinetics  of 
surface  reaction  models  which  serve  as  idealization  of  catalytic  reactions.  We  determined 
the  kinetics  of  the  saturation  process,  whereby  the  surface  becomes  filled  by  a  single  species 
and  no  further  reaction  is  possible.  We  also  examined  the  role  of  interactions  between 
adsorbates  on  the  reaction  kinetics.  Second,  we  developed  a  comprehensive  theoretical 
framework  to  describe  the  spatial  organization  of  reactants  in  two-species  annihilation, 
A  +  B  — *  0,  and  in  related  reactions.  Additionally,  we  investigated  the  reaction  dynamics 
at  an  interface  which  separates  a  population  of  A’s  and  B's.  We  determined  geometrical 
features  of  this  interface  and  the  time  dependence  of  the  reaction  rate. 

These  topics  are  treated  in  more  detail  in  the  following  report.  In  Sec.  2,  we  describe 
the  scaling  behavior  and  non-equilibrium  phase  transitions  of  surface  reaction  processes, 
and  in  Sec.  3,  we  discuss  the  wide  variety  of  interesting  features  associated  with  the  spatial 
organization  of  reactants  in  two-species  annihilation. 

2.  Kinetics  of  Catalytic  Surface  Reaction  Models 
(a)  Kinetics  of  saturation  in  the  monomer-monomer  model 

In  the  monomer- monomer  model, J  ~3  two  reactive  species,  A  and  B,  adsorb  irreversibly 
onto  single  sites  of  a  catalytic  substrate  and  nearest-neighbor  AB  pairs  can  bond  to  form 
a  reaction  product  which  desorbs  from  the  substrate.  This  bimolecular  surface  reaction  of 
the  two  adsorbed  monomeric  species  can  be  represented  by  the  following  steps: 

A  +  S  A. 

B  +  S-^B.  (1) 

A,  +  B,  -iu  ABT+2S. 

Here  A  and  B  denote  the  two  species  in  the  gas  phase,  S  denotes  a  free  surface  site,  and 
the  subscript  s  denotes  the  adsorbed  species.  If  the  adsorption  rates,  kA  and  k^,  are 
much  greater  than  the  surface  reaction  rate,  kr,  the  process  is  reaction-limited,  while  the 
process  is  adsorption-limited  in  the  opposite  case.  In  the  adsorption  step,  an  A  is  chosen 
with  probability  p  =  k^/(k^  +  &b)i  or  a  B  is  chosen  with  probability  q  =  1  —  p,  and  an 
attempt  is  made  to  adsorb  the  chosen  molecule  onto  an  empty  lattice  site.  If  adsorption 
occurs,  then  if  A,B,  nearest-neighbor  pairs  are  created,  one  such  pair  bonds  to  form  an 
AB  molecule  which  desorbs  from  the  surface  at  a  rate  governed  by  kr.  We  first  consider 
this  reaction  in  the  irreversible  limit  where  there  is  no  desorption  of  unreacted  A’s  or  B’s. 

While  this  idealized  model  does  not  precisely  correspond  to  an  actual  catalytic  reac¬ 
tion,  the  model  incorporates  the  basic  steps  of  adsorption  and  surface  reaction  that  are 
characteristic  of  most  heterogeneous  catalysis  reactions.4  Additionally,  the  simplicity  of 
the  monomer-monomer  model  renders  it  amenable  to  detailed  theoretical  analysis.  Thus 
the  investigation  of  the  monomer-monomer  model  provides  useful  insights  into  more  re¬ 
alistic  catalytic  processes.  We  first  investigated  the  kinetics  of  the  monomer-monomer 
model  in  a  mean-field  approximation,3  in  which  the  substrate  is  considered  to  be  an  N-s\te 
complete  graph.  We  calculated  the  probability  density  that  the  surface  has  a  given  density 
difference  by  analyzing  the  Fokker-Planck  equation5  that  arises  from  continuum  limit  of 
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the  underlying  master  equation  for  the  probability  distribution  of  the  particle  numbers, 
P(nA,nB)~ 

To  write  this  equation,  consider  the  reaction- limited  process  on  an  N- site  complete 
graph  in  which  each  pair  of  sites  is  connected.  In  the  reaction- limited  case,  when  an  AB 
pair  reacts,  it  is  replaced  by  AA,  BB ,  AB,  or  BA,  with  respectively  probabilities  p2,  q2 , 
pq,  and  pq,  respectively.  Thus  P(nA,nB )  evolves  according  to  a  stochastic  process  with 
the  corresponding  hopping  probabilities 
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W{riA,nB  -*  nA  ±  l,nB  T  1)  =  2  {^2  (-—■)  ,  (2) 

with  the  rate  for  the  null  process  ^(nx.ng  — »  n^.ne)  =  1  —  W(riA,nB 

W(riA,nB  — ♦  n*  — 1,*»b  +  1).  FVom  the  master  equation  for  P(n/», ns),  the  corresponding 

Fokker-Planck  equation  for  P(x )  in  the  continuum  limit,  when  p  =  q,  is 


dP(x,t) 

dt 


1  &1 

2 N  dx7 


(x(l-x)P(x,t)). 


(3) 


Here  x  denotes  the  concentration  of  the  A’s,  x  —  nA/N.  The  state-dependent  diffusion 
coefficient  D(x)  =  i(l  —  x)  reflects  the  fact  that  the  probability  of  a  reaction  event  is 
proportional  to  the  concentration  of  AB  pairs,  r(l  —  x).  Thus  the  evolution  of  the  surface 
concentration  is  analogous  to  diffusion  in  a  medium  that  is  increasingly  “sticky”  near  the 
extremities  of  a  finite  absorbing  interval,  x  -*  ±1. 

The  full  solution  to  the  probability  distribution  is  an  eigenfunction  expansion  over 
Gegenbauer  polynomials  of  order  3/2  (while  in  the  the  adsorption-limited  case,  the  expan¬ 
sion  is  over  integer-order  Bessel  functions).3  From  these  solutions,  we  determined  that 
discrete  fluctuations  ultimately  cause  the  surface  to  become  “saturated”,  i.  e.,  completely 
covered  by  only  one  species,  so  that  the  reaction  ultimately  stops.  As  the  reaction  pro¬ 
ceeds,  the  probability  that  the  catalyst  has  not  saturated  by  time  t  decays  exponentially  in 
time,  and  the  mean  time  until  saturation,  (t),  is  proportional  to  the  number  of  reactive  sur¬ 
face  sites.  Numerical  simulations  indicate  that  spatial  dimension  d  =  2  is  a  critical  value, 
above  which  the  kinetics  of  the  surface  reaction  is  accurately  described  by  the  mean-field 
behavior.  In  one  dimension,  it  is  also  obtained  an  exact  solution  in  the  rerction-controlled 
limit,6  in  which  ( t )  is  proportional  to  the  square  of  the  length  of  the  system;  this  agrees 
with  our  earlier  numerical  studies. 

When  desorption  of  reactants  is  also  allowed  to  occur  in  the  mor  jmer- monomer  model, 
a  bistability  transition  occurs  as  a  function  of  the  desorption  and  reaction  rates.1  In 
the  master  equation  for  the  evolution  of  the  concentration  difference  on  the  complete 
graph,  desorption  leads  to  both  a  diffusion  and  a  drift  term  which  drives  the  system 
toward  equal  concentration  of  A’s  and  2?’s.  This  term  counteracts  the  influence  of  the 
inhomogeneous  diffusive  term  due  to  the  surface  reaction  which  tends  to  drive  the  system 
toward  saturation.  The  corresponding  Fokker-Planck  equation  is7 


5£Sri  -  t  h  a2x  - i)PM) + «*•*» + W1  -  z)P(z't)]  ■ 

(4) 
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The  stationary  solution  to  this  equation  is  P(y )  oc  (1  —  (y/yo)2)-&,  where  yo  and  a  depend 
on  the  desorption  rate,  and  y  =  (2x  —  1)  is  the  density  difference  between  A’s  and  B's.  For 
large  desorption  rates,  P(y)  is  peaked  at  y  =  0,  (monostable),  while  for  small  rates,  P(y) 
is  peaked  at  y  =  ±1  (bistable).  Our  prediction7  of  a  transition  between  monostability  and 
bistability  has  been  observed  in  numerical  simulations.1 

(b)  Kinetics  of  interacting  catalysts  models 

We  have  investigated  the  kinetics  of  a  monomer- monomer  catalysis  model  in  which 
there  is  a  variable  excluded  volume  interaction  among  particles  of  the  same  species  (the  A’s 
for  example).8,9  For  no  interaction,  this  system  reduces  to  the  monomer- monomer  model 
for  which  there  is  a  first-order  transition  between  A -  and  J3- saturated  phases  as  a  function 
of  the  relative  adsorption  rates  of  the  two  species.  For  infinitely  strong  excluded-volume 
interactions,  one  can  view  the  process  as  involving  a  large  A  particle  whose  diameter  is 
larger  than  a  lattice  spacing  (the  “dollar”),  so  that  two  dollar  particles  cannot  be  nearest- 
neighbors,  while  the  smaller  (“dime”)  particles  have  a  sufficiently  small  radius  that  they 
do  not  exhibit  a  nearest-neighbor  exclusion.  This  geometric  constraint  on  the  dollar  par¬ 
ticles  has  far-reaching  effects  on  the  kinetics.  Instead  of  a  first-order  transition  between 
dollar-  and  dime-saturated  phases,  we  find  that  there  is  now  a  second-order  kinetic  phase 
transition  between  a  dime-saturated  phase  and  a  reactive  steady  state  as  a  function  of  the 
relative  deposition  rates  of  the  two  species.  Through  numerical  simulations,  we  concluded 
that  this  transition  is  in  the  universality  class  of  Reggeon  Field  theory,10  a  class  which 
encompasses  nearly  all  single  component  reaction  schemes  with  a  unique  absorbing  state. 

We  also  considered  the  general  case  where  the  magnitude  of  the  excluded  volume  self- 
interaction  is  variable.9  According  to  such  an  interaction,  the  probability  that  a  particle 
of  the  interacting  species  can  adsorb  onto  an  empty  site  which  is  adjacent  to  a  site  that  is 
already  occupied  by  a  member  of  the  same  species  is  reduced  by  an  amount  which  depends 
on  the  interaction  strength.  As  a  function  of  this  interaction,  the  first-order  transition  line 
for  weak  interaction  terminates  at  a  tricritical  point  where  two  second-order  transitions 
meet.  These  transitions,  which  also  appear  to  be  in  the  Reggeon  Field  Theory  universality 
class,  separate  the  dollar-saturated,  reactive,  and  dime-saturated  phases.  To  determine 
the  location  of  the  transition  and  the  corresponding  scaling  exponents,  we  have  performed 
Monte  Carlo  simulations,  series  expansions,  and  cluster  mean-field  approximations. 

From  Monte  Carlo  simulations  of  the  steady-state  and  cluster  mean-field  approxima¬ 
tions,  we  find  that  the  concentration  of  vacancies  and  dollars  both  vanish  as  the  dime 
saturated  phase  is  approached  as  c„  ~  (p  —  Pc)6*  and  cp  ~  (p  —  pc)0D ,  respectively. 
Here  p  denotes  the  relative  input  rate  of  dollars,  and  pc  is  the  value  of  this  parameter 
at  the  phase  transition.  We  thereby  estimate  the  order- parameter  exponents  for  vacan¬ 
cies  and  dollar  particles  to  be  f}v  S  0d  *=  0.28  ±  0.05.  From  subcritical  series  expan¬ 
sions  and  time-dependent  Monte  Carlo  simulations,  we  also  estimate  the  dynamical  ex¬ 
ponents  for  the  survival  probability  and  the  order  parameters  to  be  6  =  0.16  ±  0.01  and 
r)v  a  T)p  =  0.315  ±  0.005.  Here  the  survival  probability  exponent  6  is  defined  in  terms  of 
the  survival  probability,  i.  the  probability  that  the  system  has  not  yet  entered  a  satu¬ 
rated  state  at  time  t.  This  quantity  is  expected  to  vary  at  the  critical  point  as  P(t)  ~  t~s. 
Similarly  the  exponents  tjv  and  yo  are  defined  in  terms  of  the  time  dependence  of  the 
number  of  vacancies  or  dollars,  respectively,  when  the  system  is  initially  prepared  with  all 
dimes  and  a  single  vacancy.  All  of  our  exponent  estimates  are  very  close  to  those  of  the 
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Reggeon  field  theory  in  one  dimension,  suggesting  that  the  two  systems  are  in  the  same 
universality  class. 

3.  Spatial  Organization  in  Two-Species  Annihilation 
(a)  Closest  particle  to  an  absorbing  trap 

We  used  a  quasi-static  approximation  and  extreme  value  statistics  to  obtain  the  typical 
value  of  the  distance  from  an  absorbing  trap  to  the  nearest  particle,  zmi„>  in  a  Brownian 
particle  system.11  The  quasi-static  approximation  is  an  extremely  powerful  yet  simple 
theoretical  tool  to  analyze  the  spatial  distributions  of  particles  in  reacting  systems  that 
are  governed  by  diffusive  dynamics.12  In  this  method,  the  diffusion  equation  is  replaced  by 
the  steady-state  Laplace  equation,  and  the  boundary  between  the  growing  depletion  layer 
due  to  trapping  and  the  static  large- distance  profile  is  accounted  for  by  a  moving  boundary 
condition.  This  leads  to  a  greatly  simplified  form  for  the  probability  distribution  function, 
but  one  which  still  retains  all  of  the  correct  asymptotic  behavior.  We  have  then  applied 
extreme  value  statistics13  to  this  probability  distribution  to  determine  the  location  of  the 
particle  which  is  closest  to  the  trap. 

In  one  dimension,  we  find  that  the  typical  value  of  the  distance  from  the  trap  to 
the  closest  particle  grows  with  time  as  imjB  ~  t1^4,  and  that  the  probability  distribution 
for  Zinin  is  &  Gaussian  times  a  power  law.  This  reproduces  the  results  from  the  exact 
solution,  which  is  much  more  complicated,  technically.14  In  two  dimensions,  the  exact 
solution  is  extremely  involved  and  physically  unilluminating,15  while  the  quasi-static  ap¬ 
proximation  provides  a  simple  and  physically-appealing  method  of  solution.  We  find  that 
p  S  Xn, jn /(trap  radius)  grows  as  >/1n7,  and  the  controlling  factor  in  the  distribution  for  p 
behaves  as  exp(— const,  x  p7  \npfy/\nt).  Finally,  in  three  and  higher  dimensions,  we  find 
that  p  — »  const,  as  t  — ♦  oo. 

(b)  Spatial  distribution  of  reactants  in  two-species  annihilation 

In  diffusion-limited  two-species  annihilation,  A  +  B  — ►  0,  particles  organize  into  con¬ 
tinuously  growing  single-species  domains  of  linear  dimension  proportional  to  t1/2,  while 
the  typical  interparticle  spacing  grows  as  t1^4,  when  d  <  4. 16  This  spatial  organization  is 
the  mechanism  responsible  for  the  well-known  slow  kinetics,  as  reactions  take  place  only 
on  domain  boundaries  rather  than  homogeneously  throughout  the  system.  Rather  unex¬ 
pectedly,  we  found  that  a  new  length  scale,10  which  grows  as  t3/,s,  is  needed  to  describe  the 
“gaps”  between  domains  in  d  —  l.17  This  interdomain  gap  size  can  be  derived  in  terms  of 
the  known  decay  of  the  density  c(t).  If  one  postulates  that  the  gaps  are  of  length  t*B  ~ 
then  there  will  typically  be  a  reaction  in  each  gap  in  a  time  interval  At  <x  £2AB.  This  leads 
to  a  change  in  the  concentration  which  is  proportional  to  the  inverse  domain  size,  since 
one  pair  typically  reacts  per  domain.  As  a  result,  ^  ~  a  Thus  in  order  that 

c(t)  is  proportional  to  f1^4,  it  is  necessary  that  £  =  3/8. 

Concomitantly,  we  found  that  the  average  density  profile  within  a  domain  is  spatially 
inhomogeneous.  Just  as  in  the  case  of  diffusing  particle  near  a  static  trap,  the  density 
of  one  species  is  depleted  as  the  boundary  that  defines  the  opposite  species  region  in 
approached.  Particles  near  the  center  of  the  domain  are  relatively  closely  spaced,  with 
a  typical  spacing  which  grows  as  t1^4,  while  particles  are  systematically  more  distantly 
spaced  as  the  edge  of  a  domain  is  approached.  Within  an  adiabatic  approximation,  which 
is  suitable  for  slowing  growing  domains,  the  shape  of  the  density  profile  within  a  single 
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domain  has  a  sinusoidal  shape,  with  the  density  vanishing  linearly  as  the  domain  edge 
is  approached.  From  this  inhomogeneous  profile,  we  found  an  unusual  behavior  for  the 
time  dependence  of  moments  of  the  distance  distribution  between  adjacent  same  species 
particles:  Af»  s  Namely, 


(  t1/4  for  n  <  1, 

Mn  ~  \  t^lnt,  for  n  =  1,  (5) 

{  t(3n-l)/*n  for  n>l 

The  logarithmic  factor  in  the  ratio  between  the  average  and  typical  spacing  between  neigh¬ 
boring  like  particles  stems  directly  from  the  linear  decrease  in  the  density  near  the  domain 
edge.  This  same  mechanism  also  leads  to  a  power-law  diverging  factor  for  the  higher  mo¬ 
ments.  As  n  —♦  oo,  the  reduced  moment  is  dominated  by  the  sparsely  populated  region 
near  the  gap  where  the  particles  are  separated  by  a  distance  which  grows  as  t3/8 . 

Even  richer,  and  unanticipated,  spatial  organization  effects  occur  when  the  diffusion 
coefficients  of  the  two  species,  DA  and  Db,  are  unequal.  For  the  extreme  case  where 
D&  =  0,  the  spatial  correlations  of  the  mobile  and  immobile  species  are  dramatically 
different.  The  probability  of  finding  a  gap  of  length  x  between  two  nearest-neighbors  of 
the  same  species  decays  exponentially  in  x  for  the  mobile  (A)  species,  and  as  i-3''2  for  the 
static  species  (B).  The  power  law  decay  arises  because  the  spatial  correlation  of  the  B's 
is  altered  only  by  the  diffusion  of  A’s  into  B  domains.  This  kinetics  can  be  described  in 
terms  of  the  first-passage  probability  for  one-dimensional  random  walks.18  Thus  B-rich 
regions  exhibit  a  very  slow  decay  of  correlations  and  this  is  ultimately  responsible  for  the 
power-law  behavior  of  the  gap  lengths.  Similarly,  the  distribution  of  domain  lengths  L  for 
the  A’s  also  decays  exponentially  in  L,  while  the  distribution  decays  as  L~1^7  for  the  B’s. 

(c)  Dynamics  at  a  reaction  interface 

In  homogeneous  two-species  annihilation,  domains  form  spontaneously  from  a  random 
initial  reactant  distribution,  and  the  kinetics  is  governed  by  the  reaction  rate  at  the  inter¬ 
faces  between  domains.  This  motivates  our  study  of  a  system  with  a  single  well-defined 
reactive  interface,  so  that  the  microscopic  aspects  of  the  reaction  can  be  probed.19  We 
therefore  investigated  steady-state  geometrical  properties  of  the  reaction  interface  in  two- 
species  annihilation,  A  +  B  — ♦  0,  when  a  flux  j  of  A  and  B  particles  is  injected  at  opposite 
extremities  of  a  finite  domain. 

The  kinetics  of  this  process  is  described  by  the  reaction-diffusion  equations 

=  DAVhA  -  kcAc b, 

f  (6) 

=  Z?flV2ca  -  kcACB , 

where  D,  is  the  diffusion  coefficient  of  species  i  and  k  is  the  reaction  constant.  For  equal 
fluxes  of  A' s  and  B’s  injected  at  opposite  ends  of  the  system,  the  resulting  steady  state  is 
described  by  setting  the  time  derivative  to  zero  and  postulating  the  boundary  conditions 

Dc'a\zxL  =  h  Dcb\s=l  =  0  and  Dc'A\t^L  =  0,  Dc'g |tas_L  =  -j.  (7) 
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These  boundary  conditions  account  for  the  constant  and  opposing  particle  fluxes  and  the 
confinement  of  reactants  within  (-£,£]. 

Within  this  description,  we  determined  the  extent  of  the  reaction  zone  and  the  spatial 
distribution  of  the  two  species  both  by  an  asymptotic  solution  to  the  above  equations  and 
from  simple  heuristic  arguments.  For  example,  for  a  large  input  flux,  the  concentration 
profile  is  nearly  linear  near  the  domain  boundaries  with  the  magnitude  of  the  slope  pro¬ 
portional  to  jjD.  Since  the  reaction  zone  is  the  region  for  which  the  concentrations  of 
both  species  are  non-negligible,  then  the  concentration  in  the  reaction  zone  should  be  of 
the  order  jw/D ,  where  w  is  the  reaction  zone  width.  Consequently,  the  number  of  an¬ 
nihilation  events  per  unit  time  is  of  order  Arc^Cfltij,  which  is  obtained  by  integrating  the 
reaction  term  over  the  reaction  zone.  This  quantity  should  equal  the  flux  of  particles  j 
entering  the  domain.  Therefore,  balancing  these  two  rates  gives 


kcACBtV 


(8) 


or  the  reaction  zone  width  w 
reaction  zone  is 


ca,cb  ~  jw/D 


Similarly,  the  typical  concentration  in  the 


1/3 


(9) 


In  the  low  flux  limit,  the  concentration  is  almost  independent  of  position  and  is  pro¬ 
portional  to  y/j.  In  the  latter  case,  the  local  reaction  rate  reaches  a  maximum  at  the 
edges  of  the  system  rather  than  at  the  midpoint.  When  the  two  species  approach  at  a 
finite  velocity,  there  exists  a  critical  velocity,  above  which  the  reactants  essentially  pass 
through  each  other.  Results  similar  to  those  in  one  dimension  were  found  in  two-  and 
three-dimensional  radial  geometries.  Finally,  we  applied  the  quasistatic  approximation  to 
our  steady-state  solution  to  recover  the  known  time  dependence  for  the  reaction  zone  width 
for  the  case  of  initially  separated  components  with  no  external  input.  Our  approaches  are 
also  applicable  to  other  stoichiometries. 


(d)  Kinetics  and  spatial  structure  in  a  competing  population  biology  model 

We  ascertained  the  role  of  spatial  inhomogeneities20  in  the  kinetics  of  competing  bio¬ 
logical  species  by  investigating  the  kinetics  in  a  prototypical  population  dynamics  model21 
which  involves  two  competing  species.  Each  diffusing  species  undergoes  logistic  growth, 
i.  e.,  a  growth  rate  proportional  to  the  density,  and  a  self-regulation,  which  is  proportional 
to  the  density  squared.  In  addition,  there  is  a  mutual  competition  where  different  species 
tend  to  annihilate  when  they  meet.  These  processes  are  embodied  by  the  reaction-diffusion 
equations. 


—  =  D^A  -  kAA(  1  -A)-  kABAB, 

Q  D 

=  Db V2A  -  k„B(  1  ~B)~  kBAAB. 


(10) 


In  the  mean-field  limit,  where  the  species  are  assumed  to  be  homogeneously  distributed, 
there  is  coexistence  for  sufficiently  weak  competition,  kAB ,  kBA  <£  1,  while  only  one  species 
survives  when  the  competition  is  strong. 
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For  diffusing  species  in  one  dimension,  we  have  found  a  qualitatively  different  be¬ 
havior.  In  the  strong  competition  limit  in  one  dimension,  the  populations  organize  into 
single-species  domains  whose  characteristic  length  grows  logarithmically  in  time.  We  have 
determined  the  detailed  process  by  which  a  small  domain  of  B's  (for  example)  is  overrun 
by  the  larger  surrounding  domains  of  A’s.  The  smaller  domain  shrinks  rather  slowly  until 
a  critical  size  is  reached,  beyond  which  the  smaller  domain  disappears  almost  immediately 
afterward.  This  qualitative  behavior  is  also  found  to  occur  in  simulation  studies  in  two 
dimensions.  In  this  case,  domains  grow  as  a  power  law  in  time,  and  the  spatial  pattern  of 
the  domains  strongly  resembles  those  observed  in  spinodal  decomposition.22 
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